Constants of Geodesic Motion in Higher-Dimensional Black-Hole Spacetimes 
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In [l|] we announced the complete integrability of geodesic motion in the general higher- 
dimensional rotating black-hole spacetimes. In the present paper we prove all the necessary steps 
leading to this conclusion. In particular, we demonstrate the independence of the constants of mo- 
tion and the fact that they Poisson commute. The relation to a different set of constants of motion 
constructed in @] is also briefly discussed. 
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I. INTRODUCTION 

Spacetimes of higher dimensions (D > 4) have become 
much studied as a result of their appearance in theories of 
unification, such as string/M theory. Of such spacetimes, 
one important class is a sequence of black-hole metrics of 
greater and greater generality in higher dimensions that 
have been discovered over the years. 

The first such higher-dimensional black-hole spacetime 
was the metric for a nonrotating black hole in D > 4 
(the generalization of the 1916 Schwarzschild metric in 
four dimensions Q), found in 1963 by Tangherlini [J]. 
Next was the metric for a rotating black hole in higher 
dimensions (the generalization of the 1963 Kerr metric 
in four dimensions ||), discovered in 1986 by Myers and 
Perry [6j] in the case with zero cosmological constant. 
Then in 1998 Hawking, Hunter, and Taylor-Robinson 
found the general D~5 version of the D = 4 rotat- 
ing black hole with a cosmological constant (often called 
the Kerr-(anti-)de Sitter metric) that had been found in 
1968 by Carter @, Q. In 2004 Gibbons, Lii, Page, and 
Pope [l(| EH discovered the general Kerr-de Sitter met- 
rics in all higher dimensions, and in 2006 Chen, Lii, and 
Pope [l2| put these into a simple form similar to Carter's 
and were able to add a NUT [13j parameter (though not 
charge) to get the general Kerr-NUT-(A)dS metrics for 
all D. 

It is important to study the properties of these higher- 
dimensional black-hole spacetimes, and one key property 
is the nature of geodesic motion in them. In [l|, |2| we ex- 
hibited D constants of geodesic motion and announced 
that they are all independent (making the geodesic mo- 
tion integrable) and that the Poisson brackets of any pair 
of them vanish (making the integrable geodesic motion 
completely integrable) . In this paper we shall prove these 
assertions. 
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After introducing the metric and its basic symmetries, 
we recapitulate our construction of constants of geodesic 
motion and show how these constants can be generated 
from a generating function. The two main proofs demon- 
strating the independence and the Poisson commutativ- 
ity of these constants follow. The canonical formalism 
used in the text is reviewed in the Appendix. We type 
tensors in boldface with components in normal letters. 
The spacetime indices are denoted by Latin letters from 
the beginning of the alphabet, a, b, c = 1, . . . , D, and we 
use the Einstein summation convention for them. For a 
rank-2 tensor B the symbol B stands for the matrix of 
its components B a b . Where it cannot lead to a confusion 
a dot indicates contraction, i.e., a b = a e b e . We assume 
automatic lowering and raising of indices using the met- 
ric. d x a stands for the coordinate vector associated with 
the coordinate x a . 



II. HIGHER-DIMENSIONAL BLACK-HOLE 
SPACETIMES 

The general Kerr-NUT-(anti-)de Sitter spacetime dis- 
covered by Chen, Lii, and Pope may, after a suitable 
Wick rotation of the radial coordinate, be written IL2I1 



n-l 



-dx 2 

M =l " -w* a-u ' ;,-u 

with n = [D/2\ and e = D - 2n. Here, Q M = X^/U^ , 



s=£[i^(E4 fe) ^) ]-^(5>«dfc), 

(i) 



u„ = n - x i ) » x » = J2 °k x i k - 2h » x » 
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The coordinates x M (fi = 1, . . . , n) correspond to radial 
and latitude directions, ipk ( k = 0, . . . , n — 1 + e) to tem- 
poral and azimuthal directions. The parameter c n is pro- 
portional to the cosmological constant, and the remaining 
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constants Ck, c and are related to the rotation param- 
eters, the mass and the NUT parameters. Hamamoto, 
Houri, Oota and Yasui [3] derived explicit formulas for 
the curvature and demonstrated that in all dimensions 
this metric obeys the Einstein equations 

R ab = (-l) n (D-l)c n9ab . (3) 

Besides the obvious spacetime isometries generated by 
the D — n Killing vectors d^ k , the spacetime possesses 
a whole set of hidden symmetries [l|, l2|, which can be 
generated from the principal (rank-2 closed) conformal 
Killing- Yano tensor discovered by Kubizhak and Frolov 
[l5j |. These hidden symmetries play the crucial role for 
the integrability of the geodesic motion. 

The metric (JXJ) can be diagonalized. Let us introduce 
the orthonormal basis one-forms 

e>1 = / dx^ , 

n-l 

n 

e 2n+1 = (-c/AW) 1 / 2 ^^)*. (4) 

fc=0 

Then we have 

D n 

g = Y,e a e a = ^ (e"e" + e A e A ) + e e 2n+1 e 2n+1 , (5) 

a—l M— 1 

and the principal conformal Killing- Yano tensor h which 
obeys the equations 

(D - l)V a h bc = g ab £ c - g ac ^ b , £ a = V c h c a , (6) 

takes the extremely simple form 

n 

h = J2 x ^ A efl ■ ( 7 ) 

In what follows we shall also use the conformal Killing 
tensor 

Q = -hh , i.e., Q ab = h ac h bd g cd , (8) 
which takes the explicit form 

n 

Q = £a^(e"e'» + eA e A), (9) 
and satisfies ~S7( a Qbc) — 9(abQc) i where 

Q a = ^ T ^(2y c Q c a + y a Q c c ). (io) 



III. CONSTANTS OF MOTION 

In [lj we have claimed that in the spacetime (O there 
are D independent constants of geodesic motion, given 
by the following quantities: (a) n — 1 observables Cj, 
j = 1, . . . , n — 1, given by traces of powers of the projec- 
tion F of the principal conformal Killing- Yano tensors h 
(cf. Eqs. (JT5J) and (JT6|) below) 

C^trtH^F 2 ^] , (11) 

(b) D — n observables pj, j = 0, . . . , D — n — 1, given by 
symmetries of the spacetime 

Pj = u- , (12) 

and (c) the square w of the (unnormalized) velocity u 

w = u ■ u = u a u a . (13) 

Moreover, these quantities commute in the sense of Pois- 
son brackets on the phase space. Here we want to eluci- 
date and prove these properties in more detail. 

We understand all mentioned quantities as observables 
(i.e., functions) on the phase space T = T*M. It is well 
known that the cotangent space T*M has a natural phase 
space structure (cf. the Appendix or [l6j]). Since we 
investigate the relativistic theory and M is a spacetime 
manifold describing also the physical temporal direction, 
the phase space r = T*M is an unphysical phase space 
which is, however, well suited for an investigation of the 
geodesic motion. Doing canonical mechanics on it allows 
us to solve the geodesic motion in an external time which 
can be identified at the end with the affine parameter of 
the studied geodesic. 

We denote the momentum variable on the cotangent 
space as u. Indeed, since the geodesic motion is gov- 
erned by the Lagrange function L = ^ u ■ u = i u a u h g ab , 
the canonical momentum can be (up to a position of the 
tensor index) identified with the (unnormalized) velocity 
u. The Hamiltonian then is 

H = = \u u = ^u a u b g ab . (14) 

We easily realize that pj defined in (fT2|) are the special 
components of momentum and that they are constants 
of motion since d^ !j are Killing vectors. The quantities 
Cj, eq. (|11|). are constants of motion because the tensor 
F, defined as 

F ab = (h ab u c + h bc u a + h ca u b ) u c , (15) 

and the square w of the velocity it, are covariantly con- 
served along the geodesic. Indeed, thanks to (J6]), for 
u c V c u a = we have u c V c F ab = 0. 

Next we express the constants of motion Cj in terms of 
the quantities related to the principal conformal Killing- 
Yano tensor h. The components of the tensor (fTS"]) can 
be rewritten as 

w^F = PhP , (16) 
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where P is the projector orthogonal to the velocity it, 
P = I — p, i.e., P b a — 8% — p% . Here we also introduced 
the projector p, 



pi = w 1 u a u b , 



(17) 



onto the direction u. Using the cyclic property of the 
trace we thus have 



Cj = {-iyw j tr[(hP) 2j ] 



(18) 



The trace of the matrix product could be viewed dia- 
grammatically as a loop formed by joined vertices (each 
with two 'legs') corresponding to matrices in the prod- 
uct. In our case the loop is formed by alternating h and 
P vertices. Substituting P = I — p we get a sum over 
all possible loops in which P is replaced either by I or 
by ~p. In the case of the identity I the correspond- 
ing vertex is effectively eliminated, and in the case of 
the one-dimensional projector p = w~ 1 uu the loop splits 
into disconnected pieces. Namely, we can use the identity 



tv(h kl ph k2 p ■ ■ ■ h k "p) = tr(h kl p) tv(h k2 p) ■ ■ ■ tr(h k <p) 



The trace in (JTSJ) thus leads to 
tr[(hP) 2j ] =tr(h 2j ) 

+ E E (-^UliH^p) 

c—l ki<---<k c i—1 



(19) 



(20) 



The sum over c is the sum over the number of 'splits' of 
the loop, the indices fcj are the 'lengths' of the split pieces, 
and the combinatorial factor N k 3 k gives the number of 
ways in which the loop of the length 2j can be split to 
c pieces of lengths ki, . . . ,k c . From the fact that the 
tensor h is antisymmetric, it follows that traces of odd 
powers of h (optionally multiplied by a projector) are 
zero. Setting ki = 2Z, and introducing the rank-2 con- 
formal Killing tensor Q from (JHJ), eq. ([2H)l thus reduces 
to 

tr[(-hPhP) j ] = tr(Q j ) 

+ ± £ (-l) c 2^...,ntr(Q^), (21) 

C=l 1 1 <-<1 C 8=1 

h+-+l a =3 

where we used N 2 f 2 i ~ ^iVj? i which follows from the 
definition of the iV's. If we define the quantities 

Wj = wti(Qip) = u ao Q a a °Q a a l ■ ■ ■ O: u" . (22) 

we finally obtain 

Cj = w j tt[(-hPhP) j ] 

= ^tr(#)+2£ J2 (- 1 ) CN Ll i wj ~ C Il w h, 



C=l li< — <J C 
h+-+l c =3 



which is eq. (17) of 0. 

Let us note that by the same argument as that leading 
to eq. (f2"0"]) . we can derive the relation for the trace of a 
power of QP, 



tv[(QPy] =tr(Q j 



E E (-^UU^Q 1 



(24) 



Comparing with eq. (|2ip , we see that we have proved the 
relation (16) of [lj, 



c=l 1 1 <---<1 C 
h+-+lc=j 



tr[(-hPhP) j ] +tr[Q J ] =2tr[(QP) j ] 



(25) 



The relation (|23p and an algorithm for computing the 
coefficients l can be derived also in a different way. 
It was mentioned in Q that the constants Cj can be 
generated from the generating function Z{fi) = log W((3): 



Z((3) = E 

3=1 



f-lV'+l R3 ™ 1 

^— ^ — —Cj = -Y ^ tvUhPfn 

2j vfi 3 ~[ 2 i 
trlog(l -y/0hP) = logdet (i-^hP) . 



(26) 



The third equality follows from the antisymmetry of h. 
Using properties of the determinant, the antisymmetry of 
h, I = P + p, and the fact that the projector p is one di- 
mensional, we can split Z({3) into two pieces (cf. eq. (2.7) 
and (2.8) of 0): 



Z({3)=logWo(/3)+logZ(P) , 
W {(3) = det(J-v^/i) = det?(l+j3Q) , 

£(/?) = det(P + (J -VWV) 

= tr((I - = tr((/ + 0)-V) 



(27) 



Equation (|23p then corresponds to the term propor- 
tional to j3 J in the power expansion of Z(j3). The first 
term of |23|) is obtained from log Wo ((3), and the sum 
over all possible splittings of the loop corresponds to 
the f3 3 term of log£(/3). Clearly, the j-th derivative 
of logS(/3) (evaluated at f3 = 0) contains the sum over 
all possible products of l-th derivatives S'^(0) which are 
proportional to Wi defined in ([22l . The factors l2 
can thus be obtained by the explicit computation of the 
derivatives of the generating function log £(/?): 



C 4 



w ] tr{Q 3 ) 
2(-w) j <P 



log(l + £(- 



k=l 



(28) 



0=0 



(23) 



Using software for algebraic manipulation we easily get 
the first five constants (sufficient for the integrability of 
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geodesic motion up through D = 13): 

C\ = w tr Q — 2w\ , 

C 2 = w 2 tr Q 2 -4ww 2 + 2w\ , 

C 3 = w 3 tr Q 3 - 6w 2 w 3 + 6w wiw 2 - 2w\ , 

C 4 = w 4 trQ 4 -8w 3 w 4 + w 2 {Awl + 8w 1 w 3 ) (29) 

— 8w w\w 2 + 2w\ , 

C 5 = w 5 tr Q 5 - 10w i w 5 + w 3 (lOw 2 w 3 + 10wiu> 4 ) 

- w 2 (10wiwl + I0wlw 3 ) + I0wwlw 2 - 2w\ . 

Taking into account the facts that the eigenvalues of 
the principal conformal Killing- Yano tensor h are given 
by the coordinates x u , cf. eq. ((TJ, respectively, that the 
eigenvalues of Q are x^, see eq. (j9]), we can write down 
an explicit form for tr Q J and Wj : 



= £**(t£ + u£) 



(30) 



(31) 



Let us also point out that on the level of the generating 
functions the relation ([25| corresponds to 



det(J - /3 hPhP) det(l+(3Q) = det 2 (l + PQP) . (32) 

It was realized in [2| that the generating function 
W{0) =expZ(/3) = Wo(/3)S(/3) actually generates an- 
other set of conserved quantities Cj by 



W(0) 



i 00 



(33) 



3=0 



which are quadratic in the velocity u. (That they 
are quadratic can be seen from the fact that Wq does 
not depend on the velocity, from eq. pip, and from 
ioE(/3) = E^ot- 1 )^/^-) The relation between W{0) 
and Z((3) implies that 



and in particular: 



j 

E 

fc=i 



c k P k 



0=0 



(34) 



Ci = 2 Cl , 
C 2 = -4wc 2 



2c\ 



C 3 = 6w c 3 — 6wcic 2 + 2c 



l • 



(35) 



C 4 = -8w 3 c 4 + 8w 2 ac 3 + Aw 2 c 2 2 - 8wc\c 2 + 2c\ , 
C 5 = 10w 4 c 5 - 10w 3 cic 4 - 10w 3 c 2 c 3 

+ 10w 2 c\c 3 + 10w 2 c 4 C2 - lQwc\c 2 + 2c\ , 

which are the inverse of the relations (3.19) of 0. 



IV. INDEPENDENCE OF CONSTANTS OF 
MOTION 

Now we can demonstrate that the quantities w, pj and 
Cj are independent at a generic point of the phase space 
r = T*M. This means that their gradients on the phase 
space are linearly independent. To prove this, it is suffi- 
cient to show that these gradients are independent in the 
vertical direction of the cotangent bundle T*M, i.e., that 
the derivatives of these quantities with respect to the mo- 
mentum u are linearly independent. To achieve this we 
will study the wedge product of the 'vertical' derivatives. 

Let us, instead of w and Cj, consider the equivalent 
set of observables (j = 1, . . . , n — 1) 



2C 4 



3 2j 
2C Q = w, 



- w 1 ~ j C J 



1 



wtrQ J + Wj 



(36) 



where dots in the first expression denote terms which 
contain w fe with k < j, cf. eqs. (|2"3"|) . ([2"5j) . 
We are interested in the quantity 1 

J = dC A • • ■ A 9C„_! A d Po A ■ • • A dpr^-i . (37) 

Due to (fT2"j) and (f3"6")) , we have dpj = d^ Jj , and 



dCj =-— (tr Q j ) u + Q J ■ u + . . . , (38) 

where dots denote linear combinations of Q k ■ u with 
& < 3\ Q l ' u represents the vector with components 
QaiQll " ' Qa l ^u ai . From the antisymmetry of the 
wedge product it follows that 

J = mA(Q-m)A---A(Q"" 1 -m)a9^ A---A9 1/jd _, i _ 1 . (39) 

(Matrix) powers Q J of the conformal Killing tensor can 
be written as 



u=l 



E 



(40) 



The second term acts on the subspace of the vectors 
spanned on d^, j . Thus, thanks to the d^ Q A ■ ■ ■ A d^^^ 
term in the wedge product, this part can be ignored 
in (|3"9"|) . Taking into account that e^e'' = d x dx^ and 
= dx^ ■ u, the substitution of (|40|) into (|39)l leads to 



J = u 1 ...u n C/0 Xl A---A9 a; „A^ A---A^ D _ l _ 1 , (41) 



The derivative df is the vector field on spacetime M with com- 
ponents df/du a , cf. the Appendix fX] (it could be written more 
explicitly as df/du). The wedge product is, strictly speaking, 
defined for (antisymmetric) forms. However, we can easily define 
the wedge product also for the vectors or lower the vector indices 
with the help of the metric to get 1-forms. 



5 



where 



U = 



E 

permutations a 
of [0, . . . , n-1] 



2(Ti 

Signer x x 



„2ct, 



= n ( 42 ) 



fi,iy— l...n 
u<jj. 



In a generic point of the phase space we have u J ^ 
and x 2 ^ ^ X 2 , (for v) and therefore J^O there, thus 
showing that the constants of motion are independent. 



V. POISSON BRACKETS 

Finally we show that the observables u>, Cj, and pj 
Poisson commute on the phase space. 

The Poisson bracket of two functions on the phase 
space r = T*M can be written as 



{A, B} = VA ■ dB - dA ■ VB 



(43) 



where VF represents an arbitrary (torsion-free) covari- 
ant derivative which ignores the dependence of F on the 
momentum u, and dB is the derivative of B with respect 
to the momentum u, cf. the Appendix. VF 1 is a 1-form 
and dF a vector field on the spacetime M, and the dot 
indicates the contraction in spacetime tensor indices. We 
use naturally the covariant derivative V generated by the 
metric on M. 

Clearly, the commutation of any observable with the 
Hamiltonian of the geodesic motion is equivalent to 
the conservation of the observable, cf. eq. (|A13j) . so we 
have 



{w, Pj } = 0, {w,Cj} = 0. 



(44) 



The Poisson bracket between observables Pj = u ■ d 1 p j 
reduces to Lie brackets of the Killing vector fields d x p j , 
which vanish because d^,. are coordinate vector fields: 



{pi,p 3 } = ■ (Vc^J u-d^.- {Vd. H 
= [8^,6^] - u = 



(45) 



The Poisson bracket of any observable with the ob- 
servable p = I ■ u linear in momentum leads to the Lie 
derivative along the vector field I, see eq. (IA12[) : 



{Ci,Pj} = £a^pi 







(46) 



Here, the Lie derivative £g^ Ci ignores the dependence 
of Ci on the momentum u, cf. the Appendix. It van- 
ishes because is a Killing vector and the definition of 
Ci respects the symmetry of the spacetime (it does not 
depend explicitly on ipj). 

Finally, it remains to evaluate the brackets {Ci,Cj}. 
To simplify the following computation, we will study 
rescaled observables 2 



Cj = (-iyvPCj = ti[(hP) 



2J1 



(47) 



The scaling H47H differs from J36H used in the previous section. 



cf. eq. p8|) . and we denote 

P = wP = wl — uu 



(48) 



Using the cyclic property of the trace, the derivative of 
Cj in the spacetime direction is 

V„C, = 2jtr[(V a / l )P(/ l P) 2 ^ 1 ] . (49) 

Here V a h is the matrix of components V a h b c of the 
covariant derivative V/i. Substituting for V a h b c from 
eq. ([6]) and using the antisymmetry of h, we obtain 



-%p-V e Cj 



t po-oubi p a l ub 2 j-2 pa2j-l 

~^ao r b 1 "' a,l r b 2 ' ' ' a 2 j-i r e 

_„ p a 2j /,f>2j-l /.fcl pdltto 

yea 2 j r b 2 j-i 023-1 ai-^b <, 

Of p a o l,bi pa! ub2j-2 pa2j-l 

— * Sao-" b t ai r b 2 ' 023-1 r e 



(50) 



For the derivative with respect of the momentum u we 
get 



±-.d*Cj = (h^PZh^PZ . . . PZ?h d v dl ) 
+ h\P%h\ 2 P% . . . P c ^hH2 3 u C2i ■ 



(51) 



Substituting ([50]) and (j5Tj) into expression (14*3"]) for 
{Ci, Cj} and using F^w 6 = 0, we find 



joij-{Ci,Cj} - 



pa2i-l h b2,- 



— Sao* hi" 1 a i ■ ■ ■ ± b 2i -i " c t"- J -d aj " 

- upt:h\ t . . . pZ;:^ 2 ^ ■ ■ ■ Pt~ lhd2 c2^ 

= . (52) 

We thus proved that the conserved quantities w, pj, 
and Cj Poisson commute with each other. Since the 
generating function Z((3) is given by power series in 
(3 with coefficients given (up to constant factors) by 
the constants Cj, then also this function (and similarly 
W(f3) = expZ(/3)) Poisson commute with w and P j, as 
well as with itself for different choices of (3: 

{Z(p 1 ),Z(fo)} = 0, {W((3 1 ),W(p 2 )} = 0. (53) 

The same is true also for the quantities Cj generated 
from W(/3) introduced in (2j. Therefore the constants 
of motion all Poisson commute (are in involution), so the 
geodesic motion is completely integrable [H, [TtJ • 



VI. SUMMARY 

We have explicitly proved the complete integrability of 
geodesic motion in the general higher-dimensional rotat- 
ing black- hole spacetimes [l2j]. The 'nontriviaP constants 
of motion are associated with the Killing tensors which 
we generated from the principal conformal Killing- Yano 
tensor. Observables Cj are quadratic in momenta and 
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correspond to rank-2 Killing tensors, whereas constants 
Cj are of higher order in momenta and correspond to 
Killing tensors of increasing rank. 

The complete integrability of the geodesic motion is re- 
lated to the issue of separability of the Hamilton- Jacobi 
equation recently accomplished by Frolov, Krtous, and 
Kubizhak [Tel ]. The relation between integrability and 
separability on a general level has been studied in the 
series of papers by Benenti and Francaviglia (see, e.g., 
[lj|) where it was demonstrated that the separability is 
possible only if all the constants of motion, correspond- 
ing to Killing vectors and rank-2 Killing tensors, Poisson 
commute. 
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APPENDIX A: COVARIANT CANONICAL 
FORMALISM ON THE COTANGENT BUNDLE 

It is textbook knowledge [l6[ that the cotangent bundle 
T*M has the natural structure of a phase space, i.e., it is 
possesses a symplectic form ft which defines the Poisson 
bracket { , }. For a base manifold M which is equipped 
with an additional geometric structure, it can be useful 
to express phase-space quantities and operations with the 
help of quantities and operations on the base manifold. 
In this Appendix we shortly review such a procedure. 3 

We call functions on the phase space r = T*M observ- 
ables, and we write F(x, u) to emphasize the dependence 
of F on the configuration variable x £ M and on the mo- 
mentum u £ T*M. 

For the base manifold M with a (torsion-free) covariant 
derivative V (in our case the spacetime manifold with the 
metric connection), it is possible to introduce the covari- 
ant derivative of an observable F(x, u) in the horizontal 
(configurational) direction of the phase space V. For any 
base manifold vector I G TM we define 



l e V e F{x,u) 



da 



F(x(a), u(a)) 



a=0 



(Al) 



3 Similarly to the main text we type tensors in bold. Optionally, we 
write here the tensors with abstract indices \'2(i |21| which help to 
indicate tensorial operations as, for example, contraction. How- 
ever, the abstract indices do not refer to any particular choice of 
coordinates. We use small latin letters for base manifold indices 
(for tensors from TM), but we do not introduce indices for the 
phase space tensors (tensors from TT). We assume implicitly 
the tensor product, i.e., ab = a Cg) b. 



where x{a) is a curve starting from x with tangent vector 
I, and u(a) is parallel transport of u along x{a). 

For a 1-form p G T*M we can also introduce the 
derivative of an observable F(x,u) in the vertical (mo- 
mentum) direction 



p e d e F(x, u) = -^-F(x, u + ap) 
da 



Q=0 



(A2) 



Thanks to the linearity of T*M, this derivative is inde- 
pendent of any additional geometrical structure. 

Derivatives l e V e F and p e d e F are derivative operators 
on r and as such they define vector fields on T, which 
we denote 4 Z e V e and p e d e . These derivatives and vector 
fields depend ultralocally on the base manifold quantities 
I and p respectively, and we can thus introduce differen- 
tials V e -F and d e F and mixed tensor quantities ^ and 
^ by 'tearing off I and p, respectively, and by 'tearing 
off' the function F. 

Clearly, V e F is the covariant derivative of the observ- 
able F (x, u) which 'ignores' the momentum u leaving it 
covariantly constant. On the other side, the derivative 
d e F 'ignores' the configuration variable x. 

For an observable F{x,u) given by a contraction of a 
spacetime tensor field f(x) with several momenta u, 



F(x,u) = f abc -(x)u a u b u 



(A3) 



the covariant derivative reduces to the standard base 
manifold covariant derivative 



V e F(x,u) = V e f abc -{x) u a u b u 



(A4) 



The momentum derivative leaves / intact 

d e F(x, u) = f ebc -(x) u b u c ■■■ + f aec -(x) u a u c ■ ■ ■ 

■ /"' :.'•!«/„»/,••• • ... . (A5) 

A general phase space observable can then be written as 
a (infinite) sum of terms of this type. 

The mixed tensor ^ is a vector field on the phase 
space (from TIT) and a 1-form on the base manifold (from 
T*M). It is actually the horizontal lift from TM to TT 
corresponding to the covariant derivative V. The mixed 
tensor |^ is a vector field on the phase space (from TT) 
and a vector field on the base manifold (from TM). It 
gives a natural identification of the cotangent fiber T*M 
with its vertical tangent space TT*M. 

The inverse symplectic form 1~2 -1 and the Poisson 
bracket can be written as 



dx du du dx 



(A6) 



We could be more explicit and write them as l e ^ and p e . 

Similarly we could write -^r- and 4^ for quantities V e F and 
d e F introduced below. However, we use such an explicit notation 
only for the mixed tensor fields and (see below) where 
the notation V e and Q e would be too brief. 
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and 



{A, B} = V e A d e B - d e A V e B 



(A7) 



They do not depend on a choice of the covariant deriva- 
tive. Indeed, if we choose another torsion-free covariant 
derivative V on M, which can be done by specifying the 
'difference' tensor T b 



V a a b = V a a b + r b c a c , 



(A8) 



the induced covariant derivative of the phase space ob- 
servables transforms as 

V a F(x, u) = V a F(x, u) + u e T e ac {x) d c F(x, u) . (A9) 

Substituting this into (|A7j) and using the symmetry 



r o 6 c = r c ° Q we find that 



{A, B} = V e A d e B - d e A V e B 



(A10) 



i.e., the Poisson bracket is independent of the choice of 
the covariant derivative. The argument for the symplec- 
tic structure is similar. 

The Poisson bracket of an observable of type (|A3j) with 
an observable p linear in momenta u, 



p(x, u) = l c (x) u c , 



(All) 



Despite the fact that we do not need them in the main 
text, let us introduce for completeness the mixed tensor 
fields D e x and Vu e dual to ^ and -§- defined by 



dx 
dx 



■D a x = 
Vu b 







d a 
du 

ga 

du 



D b x 



. 



(A14) 



Here the dot '•' indicates the contraction of the phase 
space tensor indices. 

D e x is a vector field on the base manifold M and a 
1-form on the phase manifold T. It is actually the differ- 
ential of the bundle projection x : T*M —> M. Vix e is a 
1-form both on the base manifold M and phase space T. 

These phase space 'forms' satisfy the completeness re- 
lation 



— ^ D x + — — Vu e = 8 
ox au 



(A15) 



with 5 being the identity tensor on TI\ The symplectic 
structure f2 can be written as 



n = D e xVu P Vu P D e 



(A16) 



leads, with help of (|A4|) and (|A5j) . to the Lie derivative 



{F,p}=l e V e F-d e F(VJ c )u c 

= (rv e /° 6 - - r b -v e i a - r 

= (£if ab -)u a ub--- =£iF . 



■vj 6 



u a u b ■ ■ ■ 
(A12) 



Here £if is the standard Lie derivative on M along the 
vector field I. The last equality then defines the gener- 
alized Lie derivative £iF of the phase space observable 
F along the base manifold vector field I which effectively 
'ignores' the dependence of F on the momentum u. It 
can be extended to general phase space observables by 
linearity. It can be also defined similarly to (|A1|) with 
u(a) = <f> a u given by a flow <j> a induced by the vector 
field I acting on it. £\F can be also viewed as the deriva- 
tive of the observable F along the vector field £i on T*M 
which is called the complete lift of the vector field I on M. 

Clearly, the Poisson bracket with the Hamiltonian fjj 
leads to the covariant derivative along the u direction: 



{F, H} = u e V e F 



(A13) 



Finally, if we choose the coordinate derivative 3 asso- 
ciated with a coordinate system x a on M, 



Bd x a = , Mx a = 



(A17) 



instead of the covariant derivative V, the relations (|A16|) . 
(|A6|) , and (|A7[) reduce to the standard relations in terms 
of the canonical coordinates x a ,Ub on T, namely 



and 



ft = dx e du e — du e dx e , 



fl 1 = d x ? d Uc - d Uc d x e 



_ OA 8B OA dB 

dx e du e du e dx e 



(A18) 



(A19) 



All coordinate vectors and 1-forms in (|A18|) live on the 
phase space T. 



[1] D.N. Page, D. Kubiznak, M. Vasudevan, and 
P. Krtous, Phys. R ev. Lett. 98, 061102 (2007), 
|arXiv:hep-th/061 1083] 

[2] P. Krtous, D. Kubiznak, D.N. P age, and V.P Frolov, J 
High Energy Phys. 02 (2007) 004, |arXiv:hep-th/0612029 



[3] K. Schwarzschild, Sitzungsber. deutsch. Akad. Wiss. 

Berlin, Kl. Math. Phys. Technik, 189 (1916). 
[4] F.R. Tangherlini, Nuovo Cimento 27, 636 (1963). 
[5] R.P. Kerr, Phys. Rev. Lett. 11, 237 (1963). 
[6] R.C. Myers and M.J. Perry, Ann. Phys. (N.Y.) 172, 304 



8 



(1986). 

[7] S.W. Hawking, C.J. Hunter, and M.M. Taylor-Rob inson, [16] 

Phys. Rev. D 59, 064005 (1999), | arXiv:hep-th/9811056| 
[8] B. Carter, Phys. Lett. 26A, 399 (1968). [17] 
[9] B. Carter, Commun. Math. Phys. 10, 280 (1968). 
[10] G.W. Gibbons, H. Lii, D. N. Page and C.N. Pop e, J. [18] 

Geom. Phys. 53, 49 (2005), |arXiv:hep-th/0404008| 
[11] G.W. Gibbons, H. Lii, D.N. Page and C.N. PopeTPh ys. [19] 

Rev. Lett. 93, 171102 (2004), [arXlv:hep-th/0409155| 
[12] W. Chen, H. Lii, and C.N. Pope, Cla ss. Quant. Gray. 23, [20] 

5323 (2006), |arX"iv:hep-th/0614125| 
[13] E. Newman, L. Tamburino, and T. Unti, J. Math. Phys. 

(N.Y.) 4, 915 (1963). [21] 
[14] N. Hamamoto, T. Houri, T. Oota, and Y. Y asui, J. Phys. 



[15] 



A40, F177 (2007) _ 
D. Kubiznak and V.P 



arXiv:hep-th/0611285 



Fl (2007), |arXiv:gr-qc/0 610144 

V.I. Arnol'd, Mathematical Methods of Classical Mechan- 
ics (Springer- Verlag, Berlin, 1978). 

V.V. Kozlov, Usp. Mat. Nauk 38:1, 3 (1983) [Russ. 
Math. Surv. 38:1, 1 (1983)]. 

V.P. Frolov, P. Krtous, and D. Kubiznak, J. High Energy 

Phys. 02, 005 (2007), [axXiv:hep-th/0611245l 

S. Benenti and M. Francaviglia, Gen. Rel. Grav. 10, 79 

(1979). 

R. Penrose and W. Rindler, Spinors and Space-Time 
(Cambridge University Press, Cambridge, England, 1984, 
1986). 

R. M. Wald, General Relativity (The University of 
Chicago Press, Chicago and London, 1984). 



Frolov, Class. Quant. Grav. 24, 



